Effects of phase transitions in devices actuated by the electro- 
magnetic vacuum force 
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Abstract. - We study the influence of the electromagnetic vacuum force on the behaviour of a 
model device based on materials, like germanium tellurides, that undergo fast and reversible metal- 
insulator transitions on passing from the crystalline to the amorphous phase. The calculations are 
performed at flnite temperature and fully accounting for the behaviour of the material dielectric 
functions. The results show that the transition can be exploited to extend the distance and 
energy ranges under which the device can be operated without undergoing stiction phenomena. 
We discuss the approximation involved in adopting the Casimir expression in simulating nano- 
and micro- devices at finite temperature. 



Introduction. — With the development of modern 
technology towards smaller structures an increasing at- 
tention has been addressed to the role of electromagnetic 
vacuum fluctuation forces (dispersion or van der Waals 
and Casimir forces) in micro- and nano-electromechanical 
systems [1-7]. These forces vary typically with the third 
or fourth inverse power of the distance between the sur- 
faces of the interacting bodies and therefore can be very 
intense in the sub-micrometer regime. When the moving 
parts of a device come to such a close distance, disper- 
sion forces affect the dynamics and the operation of the 
device and, under certain circumstances, they can cause 
adhesion or stiction between the surfaces, thus limiting 
the device lifetime. Devices actuated by Casimir force 
have been designed and demonstrated [8-10]. Since these 
forces depend upon the reflectivity, the size and the geom- 
etry of the interacting parts [11,12], several possibilities 
can be explored to devise structures that exploit such de- 
pendence in order to improve the performance of existing 
devices or to design new types of devices [13]. 
In this paper we investigate the possibility of using ma- 
terials, that undergo fast and reversible phase transitions, 
to extend the travel range of a device. Thermally induced 
phase changes between the amorphous and the polycrys- 
talline state in a thin film have been observed since a 
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Fig. 1: (a) Three layers model of equation (1); (b) Simple 
representation of the dispersion forces based device. 



long time in tellurides containing Ge and Sb (GST ma- 
terials) [14-21]. The transition is accompanied by signifi- 
cant changes in optical and transport properties, a feature 
that is exploited in optical data storage [22] and could be 
useful in new nanoscale memories [23,24]. Previous the- 
oretical studies have shown that metal insulator transi- 
tions may cause significant changes in the Casimir force 
behaviour [25,26]. In this paper we propose to use this 
transition to modify the force between the components of 
a simple device. We show that the locations of the ex- 
trema in the total energy potential curve of the device are 
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displaced by the transition in such a way that the range 
of distances and the energy interval over which the device 
can be operated can be significantly modified. 

The model device. — We consider an ideal actua- 
tor [1, 13] consisting of two parallel plates separated by a 
gap, with one plate fixed on a substrate and the other sus- 
pended by an elastic restoring force F — —Kxq6, where xq 
is the unactuated distance and K is the spring constant. 
The stationary plate has a flat surface at S — 1, while 6 = 
denotes the equilibrium position of the movable plate in 
the absence of dispersion forces and corresponds to the 
unstretched state of the spring. The attractive interaction 
between the plates of Fig. 1(a) is provided by the elec- 
tromagnetic vacuum fluctuation force, whose expression 
in terms of the temperature T, the dielectric functions ci 
and £2 and the inter-plate distance d = (1 — ^)a:o, is given 
by the Lifshitz formula [27,28]: 
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here A is the plate surface, (3 = l/fcsT, ks is the Boltz- 
mann constant, ri„ = 2Tm/hp is the Matsubara frequency 
corresponding to the n-th thermal fluctuation mode, the 
prime on the summation indicates that the n = term is 
given half weight. We assume the plate size to be much 
larger than the separation distance. Qtm and Qte refer 
to transverse magnetic (TM) and transverse electric (TE) 
modes respectively and are given by: 
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and the dielectric functions are evaluated at the frequency 
i^n by means of the London transform: 
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where S[e(y)] is the imaginary part of the dielectric func- 
tion. In the calculation for the metallic phase we follow 
the prescription of ref. [29], i.e. we put the n = contri- 
bution of the TE modes equal zero. The behaviour of the 
force as a function of d is non linear: at short distances it 
can be reproduced by a third inverse power of d (van der 
Waals regime), but a single inverse power term does not 
reproduce its behaviour for d values larger than a few tens 
of nanometers. As pointed out by several authors [30,31], 
mainly in connection with electrostatic actuators, there 




Fig. 2: Imaginary parts of amorphous and crystalline GeTe 
dielectric functions and their London transforms, adapted from 
ref. [36] as discussed in the text. 



is an intrinsic instability in a device of this sort, which 
prevents the plates to be stably positioned over a large 
distance. The system is bistable, having a total potential 
energy with a local and an absolute minimum separated 
by a barrier. This limits the range of motion of the de- 
vice and determines stiction at short distances between the 
plates [1,2,4-6]. We adopt the configuration illustrated in 
Fig. 1(b) where both plates are made of telluride films 
deposited on proper substrates. For simplicity we assume 
that the thickness of the deposited films is larger than the 
distance between the plates (if the film size is comparable 
to the distance, the substrate is expected to affect the in- 
tensity of the force). This assumption can be removed by 
extending the expression of the Lifshitz force to deal with 
a five layer system [13,32]. 

GeTe dielectric properties and the force. To il- 
lustrate our model device, we focus our attention on GeTe, 
whose optical properties have been experimentally deter- 
mined in both phases. It undergoes rapid transitions be- 
tween polycrystalline and amorphous states under either 
optical or electrical excitations. According to the band 
theory in a perfect crystal rocksalt structure GeTe should 
exhibit a semiconducting behaviour. However the sta- 
ble polycrystalline phases are characterized by large va- 



cancy concentrations (typically 10 



') and local dis- 



tortions [33,34]. Compared to other defects, that may 
be present, germanium vacancies have the lower forma- 
tion energy [35] . GeTe thin films show a p-type electrical 
conductivity with a high carrier density of up to 5 x 10^*^ 
cm^^. This density is consistent with the concentration 
of Ge vacant sites, i.e. the metallic conductivity is the 
consequence of the holes that are formed per single va- 
cant site. In the optical properties direct evidence of the 
metallic behaviour is provided by the presence of Drude- 
like intraband absorption in crystalline phase [36] (see Fig. 
2). On the other hand the amorphous phase has typical 
semiconductor properties, as illustrated by the behaviour 
of the imaginary part of the dielectric function, reported 
in the same figure. Conventional semiconductors, such as 
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Fig. 3: Relative difference between -Fee, Fca and Faa- 

Si or GaAs, do not show such large differences in the opti- 
cal properties as a consequence of the amorphous-crystal 
transition. 

In Fig. 2 we display the calculated London transform for 
both phases. To obtain this curve we have reproduced 
the experimental data available in the range of frequen- 
cies 0.1-^5 eV with a Drude-Lorentz model including both 
intra- and inter-band transitions in the metallic phase, and 
pure inter-band transitions in the amorphous phase. The 
values of the carrier density, the carrier effective mass, 
the relaxation time and the plasma frequency agree with 
those derived from various experimental studies [15-17]. 
It is seen that t{iu)) turns out to be very different in the 
two phases in almost all the range of frequencies. As a 
consequence of the intraband absorption, that is present 
in the metallic phase only, the curve for the polycrystalline 
phase for small w takes much higher values than the one 
for the amorphous phase. Such differences over all the 
integration interval of equation (1) are expected to deter- 
mine significant changes in the force. We can use equation 
(1) to calculate the force between two amorphous coated 
plates Faa, two crystal coated plates Fee-, and the mixed 
configuration Fca- In Fig- 3 we report, as a function of 
the separation distance, the relative variation of the dis- 
persion force between the plates with respect to the con- 
figuration with both plates in the crystalline phase, when 
one or both plates undergo the transition from the crys- 
talline to the amorphous phase. The curves have a similar 
behaviour with two ranges, at small and large distances, 
where the relative difference is constant, and a maximum 
in the range between 50 nm and 5 [im. The constant 
values are a consequence of the fact that the forces have 
identical behaviour as a function of the plate separation. 
At short distances the force can be approximated by [37]: 
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which gives a relative change of approximately 28%, that 
does not depend upon the distance and the temperature. 
Notice that for the n = term of the Matsubara sum, the 



crystalline dielectric function diverges and: 
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where C,(x) is the Riemann zeta function. On the other 
hand, at large distances within the present treatment of 
the temperature dependence of the Casimir force, the force 
between polycrystalline plates is simply given by [29,38]^: 
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In the amorphous case the force for large d cannot be 
written in a closed form. However it can be very well 
approximated by the n = term of the Matsubara sum 
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where Lin[x] is the polylogaritmic function of n-th order 
in the argument x. From these expressions one obtains a 
relative change of the order of 36%, in agreement with the 
exact numerical results in Fig. 2. 

Tailoring the device performance. Device insta- 
bility is usually discussed in terms of a dimensionless pa- 
rameter, like the ratio between the Casimir energy at the 
equilibrium distance and the elastic energy [?,1,13]. This 
can be done conveniently when the force has a simple in- 
verse power dependence with the separation distance. The 
analysis of the bifurcation diagrams as a function of the 
parameter allows one to draw conclusions about changes 
in the critical separation with the material dielectric prop- 
erties. In our case the force as a function of the separation 
distance cannot be reproduced by a simple power law, so 
that we cannot adopt a single dimensionless parameter to 
study the device instability. 

The behaviour of the device depends upon the spring con- 
stant k, the area A of the plates and the unactuated dis- 
tance xq. The ratio E = A/k gives a measure of the rela- 
tive importance of the vacuum and the elastic force, since 
they increase linearly with A and k respectively. Accord- 
ing to the values attributed to these parameters one ob- 
tains different intervals of distances and energies of device 
operation at fixed xq. To show the effects of the phase 
transition we give in Fig. 4 plots of the potential energy 
of the device for fixed k and A at xq values of the or- 
der of hundred nanometers. This choice corresponds to 
distances where the relative difference of the force varies 
significantly (see Fig. 3). As expected, on passing from 
the CC to the AA configuration the weakening of the vac- 
uum force shifts the stable equilibrium position towards 
smaller S values and the maximum at larger S. The dis- 
placement is more pronounced for low S i.e. soft spring 



^Notice that this behaviour is different from the expression of the 
Casimir force at T = 0° _fC that for a perfect metal gives Fqq = 

— 240^^4 ■ This formula has been frequently used in previous studies 
of the role of electromagnetic forces in stiction phenomena. 
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Fig. 4: Potential profiles for CC configuration (dashed line), 
CA configuration (continuous line) and AA configuration (dot- 
dashed line). In the (a) plot xo = 100 nm, fc = 0.1 N/m and 
A = 2- 10"^° letters A and B refer to Fig. 5(a); in the (b) 
plot xo = 200 nm, k = 0.5 N/m and A = 10"** 



and small area (Fig. 4(a)), where it can be of the order 
of 80% and more, than in the case of large S plotted in 
Fig. 4(b), where the increase of the travel range is of the 
order of 10%. The range of energies for which the device 
can be operated is modified in a similar way. These curves 
provide examples of the possibility of tuning in reversible 
way the device properties using the metal-insulator phase 
transition. Fig. 5 displays bifurcation diagrams of the de- 
vice as a function of S for xq equal to 100 nm (Fig. 5(a)) 
and 5 ^J^m (Fig. 5(b)) [13,39,40] for the CC, CA and 
AA configurations. The curves have similar shapes, with 
a lower branch before the fold that corresponds to stable 
solutions while the upper branch indicates the unstable 
states. The 5 values which allow to operate the device at a 
given value S* are those which lie within the intersections 
of the line 5 = 2* with the upper and the lower branch of 
the curves (points A and B in Fig. 5(a)). At the S value 
corresponding to the curve fold there is only one intersec- 
tion, which gives the critical separation Sq. As expected, 
the crystalline-amorphous transition changes significantly 
the range of S parameters over which the device can work. 
Notice that 5q varies between 0.21 and 0.22. The critical 
value is expected to change when the functional depen- 
dence of the vacuum force with the plate separation is 
modified: for the dependence typical of the Casimir 
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Fig. 5: Bifurcation diagrams for CC configuration (dashed 
line), CA configuration (continuous line) and AA configuration 
(dot-dashed line). In the (a) plot = 100 nm, the dotted line 
represents the H value used in Fig. 4(a). In the (b) plot a;o = 5 
/im, the dotted line represents the force calculated using the 
T = 0° K Casimir's force. 



force i5o = 0.20, while for the behaviour, which is char- 
acteristic of the van der Waals interaction and of the large 
distance force at T > 0° K, 5q = 0.25 [1, 13]. In the case 
illustrated in Fig. 5(a), i.e. 100 nm distance, the critical 
value deviates significantly from the third inverse power 
behaviour. The curves of the 5 fj,m case show a similar 
behaviour. For the sake of comparison we have reported 
in Fig. 5(b) the curve corresponding to a pure Casimir in- 
teraction, which has been often adopted in similar device 
analysis. This force is appropriate to describe the case of 
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to avoid inconsistencies related to the Drude description 
of the intraband absorption [43]. 
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Fig. 6: Jo as a function of xq. Dashed line corresponds to 
the T = 0° K Casimir's force, dot-dashed line is obtained in 
the small and large d limits of equations (5) and (7). Dots 
represent the exact results for Fee (circles), Faa (triangles) 
and FcA (stars). 



metallic plate at large distances for T ^ 0° K. Compari- 
son with the results for our CC configuration shows that 
the ranges of S and i5 values are considerably underesti- 
mated by this interaction. Notice that the vacuum force 
at such distance cannot yet be described by a third inverse 
power law, as expected from equation (7) at T > 0° K. 
This is most clearly seen by the critical Sq value which is 
smaller than 0.25. 

To better illustrate how the separation distance depen- 
dence of the vacuum force can affect the device perfor- 
mance, we plot in Fig. 6 the critical value as a function of 
the unactuated distance. It is seen that So varies between 
0.25 at very short (less than 10 nm) and very large (more 
than 10 fim) distances and a minimum slightly larger than 
0.20. The minimum occurs at different distances in the 
three configurations and it is found at smaller xq for the 
AA case. It is clear from this curve that simulating device 
behaviour on the basis of simple inverse power law forces 
does not allow to get realistic results. 

Conclusions. — We have shown how the metal- 
insulator transition in germanium tellurides can be ex- 
ploited in a simple device to modify its performance and 
the conditions under which it can be used. The choice 
of the appropriate values for the parameters entering into 
our model may depend upon a number of factor, like the 
film quality [41], the possibility of controlling the kinetics 
of the phase transition, the mechanical properties of the 
components, the surface roughness [42] etc. which have 
to be clarified in order to plan the realization of a spe- 
cific device. On the theoretical side two aspects have to 
be further investigated. The first concerns the role of the 
film size in determining the change in the vacuum force: 
this can be carried out along the lines of previous stud- 
ies [13, 32] with the purpose of determining the minimum 
thickness that allows to detect significant vacuum force 
variations. The second has to do with the effects of the 
finite plate area, which in the metal case may be needed 
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